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METHOD AND SYSTEM OF GENERATING AN MRS SPECTRUM FROM 

MULTIPLE RECEIVER DATA 

BACKGROUND OF THE INVENTION 

5 The present invention relates generally to magnetic resonance spectroscopy 

(MRS) and, more particularly, to a system and method for multiple-receiver proton 
spectroscopy such that a single absorption spectrum is generated as a combination of data 
received from multiple receiver coils, 

10 When a substance such as human tissue is subjected to a uniform magnetic field 

(polarizing field Bo ), the individual magnetic moments of the spins in the tissue attempt 
to align with this polarizing field, but precess about it in random order at their 
characteristic Larmor frequency. If the substance, or tissue, is subjected to a magnetic 
field (excitation field Bi) which is in the x-y plane and which is near the Larmor 

15 frequency, the net aligned moment, or "longitudinal magnetization", M z , may be rotated, 
or "tipped", into the x-y plane to produce a net transverse magnetic moment M t . A signal 
is emitted by the excited spins after the excitation signal Bi is terminated and this signal 
may be received and processed to form an image or a magnetic resonance spectroscopy. 

20 When utilizing these signals to produce images, magnetic field gradients (G*, G y , 

and G z ) are employed. Typically, the region to be imaged is scanned by a sequence of 
measurement cycles in which these gradients vary according to the particular localization 
method being used. The resulting set of received NMR signals is digitized and processed 
to reconstruct the image using one of many well known reconstruction techniques. 

25 

Magnetic Resonance Spectroscopy (MRS) may be used in vivo for the 
determination of individual chemical compounds located within a volume of interest. 
The underlying principle of MRS is that atomic nuclei are surrounded by a cloud of 
electrons which slightly shield the nucleus from any external magnetic field. As the 



1 



Docket No. GEMS808L151 



structure of the electron cloud is specific to an individual molecule or compound, the 
magnitude of this screening effect is then also a characteristic of the chemical 
environment of individual nuclei. Since the resonant frequency of the nuclei is 
proportional to the magnetic field it experiences, the resonant frequency can be 
5 determined not only by the external applied field, but also by the small field shift 
generated by the electron cloud. Detection of this chemical shift, which is usually 
expressed as "parts per million" (PPM) of the main frequency, requires high levels of 
homogeneity of the main magnetic field Bo. 

10 Typically, MR proton spectroscopy is used to generate a one-dimensional (ID) 

frequency spectrum representing the presence of certain chemical bonds in the region of 
interest. In medical diagnosis and treatment, MRS provides a non-invasive means of 
identifying and quantifying metabolites from a region of interest, often the human brain. 
By finding the relative spectral amplitudes resulting from frequency components of 

15 different molecules, medical professionals can identify metabolites indicative of diseases, 
disorders, and other pathologies such as Alzheimer's disease, cancer, stroke, and the like. 
In this context, two nuclei are typically of particular interest, H and P . Phosphorus 3 1 
MRS is directed to the detection of compounds involved in energy metabolism relating to 
membrane synthesis and degradation. Metabolites of particular interest in proton MRS 

20 studies include glutamate (Glu), glutamine (Gin), choline (Cho), creatine (Cre), N- 
acetylaspartate (NAA), and the inositols (ml and si). 

Recently, increased sampling resolutions and the ability to simultaneously receive 
data from multiple-channel receiver coils has lead to greater accuracy in MR data 
25 acquisition. Accordingly, radiologists and other medical personnel are able to review 
separate MRS results from each channel for data acquired during a single scan. 
Additionally, improved data acquisition capabilities have created opportunities for 
implementing new signal processing algorithms to improve the sensitivity and accuracy 
of MRS procedures. Specifically, nonparametric techniques, such as Capon analysis and 
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Amplitude and Phase Estimation of a Sinusoid (APES) analysis can provide 2D results 
for MRS experiments that show both frequency and damping information. These 
techniques provide improved sensitivity and differentiation between metabolites. 

5 However, whether implementing ID or 2D proton spectroscopy, interpreting the 

results of multiple channels can be a tedious and time consuming process because a 
medical professional must inspect the results from each channel independently. That is, 
since a separate result for each channel is generated, when reviewing the results, the 
medical professional must independently interpret the results from each channel. This 

10 process of reviewing multiple sets of results and changing therebetween creates an 
arduous review process that is susceptible to erroneous human estimations and decreases 
patient throughput. That is, since the medical professional reviewing the results must 
repeatedly switch between the results from each independent channel, the medical 
professional is required to estimate an overall review of data acquired from the various 

15 independent channels. For example in an 8 channel head coil, a medical professional 
must individually review data from 8 different receiver channels to make a diagnosis. 
Relying on human estimation not only lengthens the review process but also complicates 
the diagnostic process. 

20 It would therefore be desirable to have a system and method capable of generating 

a single set of results from a multiple receiver MR proton spectroscopy scan. 
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BRIEF DESCRIPTION OF THE INVENTION 

The present invention provides a system and method of generating a single 
compilation of results from an MR proton spectroscopy scan that overcomes the 
aforementioned drawbacks. Specifically, the present invention provides a system and 
5 method of processing MR signals individually and then combining the multiple sets of 
MRS results to generate and display a single compilation of MRS results. As such, a 
radiologist or other medical personnel may make a diagnosis from a single composite 
spectrum. 

10 In accordance with one aspect of the invention, a method of multi-channel MR 

spectroscopy (MRS) is disclosed that includes simultaneously acquiring MR signals from 
multiple coils and processing the MR signals individually to generate multiple sets of 
MRS results. The method also includes combining the multiple sets of MRS results to 
form a single superset of MRS results and generating and displaying a proton MRS 

15 absorption spectrum from the single superset of MRS results for clinical inspection. 

In accordance with another aspect of the invention, a magnetic resonance imaging 
(MRI) apparatus includes a plurality of gradient coils positioned about a bore of a magnet 
to impress a polarizing magnetic field and is used to generate a set of MRS results. An 

20 RF transceiver system and an RF switch are controlled by a pulse module to transmit and 
receive RF signals to and from an RF coil assembly to acquire MR images. The MRI 
apparatus also includes a computer programmed to acquire MR signals from multiple 
receiver channels and process the MR signals from each channel independently. The 
computer is further programmed to process MR signals from each coil into a single set of 

25 MRS results for each coil, then combine these results into a single composite set of MRS 
results and display a proton MRS absorption spectrum from the single composite set of 
MRS results. 
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In accordance with another aspect of the invention, the invention is embodied in a 
computer program stored on a computer readable storage medium and having instructions 
which, when executed by a computer, cause the computer to process MRS results from 
the multiple coils independently. The computer combines the processed MRS results into 
5 a single composite set of MRS results and displays the set as a single MRS absorption 
spectrum. 

Various other features, objects and advantages of the present invention will be 
made apparent from the following detailed description and the drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The drawings illustrate one preferred embodiment presently contemplated for 
carrying out the invention. 
In the drawings: 

5 Fig. 1 is a schematic block diagram of an MR imaging system for use with the 

present invention. 

Fig. 2 is a flow chart setting forth the steps of compiling a single absorption 
spectrum from a multiple receiver MR proton spectroscopy scan in accordance with one 
embodiment of the present invention. 
10 Fig. 3 is an example of a single ID absorption spectrum compiled from a multiple 

receiver MR proton spectroscopy scan in accordance with one embodiment of the present 
invention. 

Fig. 4 is an example of a single 2D absorption spectrum compiled from a multiple 
receiver MR proton spectroscopy scan in accordance with one embodiment of the present 
15 invention. 



6 



Docket No. GEMS808 1.151 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Referring to Fig. 1, the major components of a preferred MRI system 10 
incorporating the present invention are shown. The operation of the system is controlled 
5 from an operator console 12 which includes a keyboard or other input device 13, a 
control panel 14, and a display screen 16. The console 12 communicates through a link 
18 with a separate computer system 20 that enables an operator to control the production 
and display of images and MRS results on the display screen 16, The computer system 
20 includes a number of modules which communicate with each other through a 

10 backplane 20a. These include an image processor module 22, a CPU module 24 and a 
memory module 26, known in the art as a frame buffer for storing image data arrays. The 
computer system 20 is linked to disk storage 28 and magneto optical disk (MOD) drive 
30 for storage of image data and programs, and communicates with a separate system 
control 32 through a high speed serial link 34. The input device 13 can include a mouse, 

15 joystick, keyboard, track ball, touch-activated screen, light wand, voice control, or any 
similar or equivalent input device, and may be used for interactive geometry prescription. 

The system control 32 includes a set of modules connected together by a 
backplane 32a. These include a CPU module 36 and a pulse generator module 38 which 

20 connects to the operator console 12 through a serial link 40. It is through link 40 that the 
system control 32 receives commands from the operator to indicate the scan sequence 
that is to be performed. The pulse generator module 38 operates the system components 
to carry out the desired scan sequence and produces data which indicates the timing, 
strength and shape of the RF pulses produced, and the timing and length of the data 

25 acquisition window. The pulse generator module 38 connects to a set of gradient 
amplifiers 42, to indicate the timing and shape of the gradient pulses that are produced 
during the scan. The pulse generator module 38 can also receive patient data from a 
physiological acquisition controller 44 that receives signals from a number of different 
sensors connected to the patient, such as ECG signals from electrodes attached to the 

30 patient. And finally, the pulse generator module 38 connects to a scan room interface 
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circuit 46 which receives signals from various sensors associated with the condition of 
the patient and the magnet system. It is also through the scan room interface circuit 46 
that a patient positioning system 48 receives commands to move the patient to the desired 
position for the scan. 

5 

The gradient waveforms produced by the pulse generator module 38 are applied 
to the gradient amplifier system 42 having G x , G y , and G z amplifiers. Each gradient 
amplifier excites a corresponding physical gradient coil in a gradient coil assembly 
generally designated 50 to produce the magnetic field gradients used for spatially 

10 encoding acquired signals. The gradient coil assembly 50 forms part of a magnet 
assembly 52 which includes a polarizing magnet 54 and a whole-body RF coil 56. A 
transceiver module 58 in the system control 32 produces pulses which are amplified by 
an RF amplifier 60 and coupled to the RF coil 56 by a transmit/receive switch 62. The 
resulting signals emitted by the excited nuclei in the patient may be sensed by the same 

15 RF coil 56 and coupled through the transmit/receive switch 62 to a preamplifier 64. The 
amplified MR signals are demodulated, filtered, and digitized in the receiver section of 
the transceiver 58. The transmit/receive switch 62 is controlled by a signal from the 
pulse generator module 38 to electrically connect the RF amplifier 60 to the coil 56 
during the transmit mode and to connect the preamplifier 64 to the coil 56 during the 

20 receive mode. The transmit/receive switch 62 can also enable a separate RF coil (for 
example, a surface coil) to be used in either the transmit or receive mode. 

The MR signals picked up by the RF coil 56 are digitized by the transceiver 
module 58 and transferred to a memory module 66 in the system control 32. A scan is 
25 complete when an array of MR raw data has been acquired in the memory module 66. 
This MR raw data is input to an array processor 68 which processes the MR raw data as 
necessary to create the desired results which may include image data or MRS results. 
The image data or MRS results are conveyed through the serial link 34 to the computer 
system 20 where they are stored in memory, such as disk storage 28. In response to 
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commands received from the operator console 12, the image data or MRS results may be 
archived in long term storage, such as on the MOD drive 30, or may be further processed 
by the image processor 22 and conveyed to the operator console 12 and presented on the 
display 16. 

5 

The present invention includes a system and method suitable for use with the 
above-referenced MR system, or any similar or equivalent system for obtaining MR data. 
Proton MRS can be implemented on the above-referenced MR system and is used in vivo 
to measure the concentration of a number of metabolites. High-field magnetic resonance 
10 scanners (1.5T or greater) are typically used to perform MRS studies involving a number 
of pathologies, including but not limited to, cancer, multiple sclerosis, Alzheimer's 
disease, and acute stroke. 

Referring now to Fig. 2, an MRS scan is conducted to simultaneously acquire data 
15 for a ID or 2D absorption spectrum plot from multiple receiver channels of an MR 
system. The MRS scan includes a reference data acquisition in which a signal is acquired 
from a region or volume of interest. The MRS scan also includes a water-suppressed data 
acquisition in which the dominant water signal has been removed via a chemical selective 
saturation, such as a known CHESS pulse sequence. The MRS data acquisition requires 
20 that a number of scan parameters be determined for each region of interest. The center 
frequency, flip angle, transmit gain, receive gain, and shimming parameters are 
determined for each MRS experiment during prescan. These parameters are optimized 
for a particular TE and TR selected for a particular scan. 

25 After sampling the acquired MR signals, the MR data can be represented as a 

mixture of R complex damped sinusoids: 

An] = X 4e y V" ff ' +>> + a [n] (Eqn. 1) 
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where A r ,6 r ,cr r , and co r are the amplitude, phase, damping factor and angular frequency 
of the V th component, respectively, and^[«] is a noise term. 

After data acquisition, the raw data is pre-processed as ID data. That is, once 
5 non- water suppressed reference data is collected 114, averaging is used to yield an 
averaged set of reference data, r[«], where nis an index representing each sampled 
complex data point. This raw data averaging is done for each coil. 

The reference data, r[n], undergoes DC mixing 116, zero phasing 118, linear 
10 phase correction 120 and phase spline smoothing 122, to determine a phase correction 
vector 124, c[n]. Additionally, the reference data, r[n], is used to determine weighting 
125 for the data of a respective coil. Simply, weighting 125 is carried out to determine 
data influence from a particular coil on a combined spectrum or image. Specifically, if it 
is determined that the data of a particular coil should be used, the reference data, r[/z], is 
15 also used to assign weight to the receiver channels. 

The weighting 125 can be determined through multiple considerations. For a 
particular coil, a maximum magnitude for the averaged reference data, can be used 
as a particular weighting factor, w x . Additionally, channels that have a significantly weak 

20 signal relative to the channel with the strongest signal are not used for the combined 
results and assigned a weight factor of zero. Specifically, in multi-channel MRS, regions 
of interest close to receive coil elements benefit from improved signal-to-noise ratios. 
Therefore, data received from regions of interest close to receive coil elements are 
appropriately weighted to take advantage of the improved signal-to-noise ratios. The 

25 signal characteristics from each coil depend on a number of factors including the 
orientation of the coil with respect to the Bo field, the proximity of the coil to the volume 
generating the signal, coil loading, coil-to-coil coupling effects, and the permeability and 
permitivity of the medium through which the radio-frequency signal travels prior to being 
received by the coil elements. Accordingly, the current invention weights data received 

10 



Docket No. GEMS8081.151 



to take advantage of the improved signal-to-noise ratios. Once weighted, the data is then 
normalized according to: 

1.0 = 2^/ (Eqn. 2); 

L 

where L refers to the number of channels used. 

5 

Water-suppressed data is also collected from each coil element and the data from 
each coil element is averaged 112 to obtain an averaged water-suppressed data set, s[n]. 
A phase correction vector 124, c[«], is then applied 127 to the reference data set, r[n], 
and the water-suppressed data set, s[n] at 126. The phase correction vector, c[n], may be 
10 defined by: 

c[n] = • g[n] ' w[n] (Eqn. 3); 

wherein w[n] is a windowing function, A 2p is defined by r*[0], and g[n] is defined by: 

g [n] = e '««"^>**M> (Eqn. 4); 

wherein o m is the largest frequency component of the non-water suppressed data set, 
15 s[n], cojp is the number of phase cycles in the non- water suppressed data set, s[n], and 
O s [n] is the phase spline smoothing factor from an unwrapped phase of: 

r[ny K ^ )n (Eqn. 5). 

20 Applying the phase correction vector, c[n] 9 removes residual water signal from the water- 
suppressed data, s[n], and also removes distortions due to eddy currents. Specifically, 
the correction vector, c[ri]> is multiplied by the reference data set, r[ri], and the water- 
suppressed data set, s[n] as: 

r C orrected[n] = c[n] - r[n] (Eqn. 6) 

25 and 

Scorrected[n] = c[n] - s[n] (Eqn. 7), 
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respectively. Water subtraction 128, in which a scaled version of the reference data set is 
subtracted from the water-suppressed data set of each receiver channel, is performed and 
a Fourier Transform is applied 130 to the phase-corrected, water-suppressed data set, 
S CO rrected[ri\, with residual water removed. Prior to computing a frequency spectrum, 
5 SresuitsM, the data is zero-padded for greater resolution. The results of the Fourier 
Transform 130, for each receiver coil, where the number of receive coils is L, is 
combined 132 by summing the computed results from the frequency spectrum Sresuits[k] , 
wherein 

Sresuits(k) = Re{FT{ssignai[n]}} (Eqn. 8) 

10 and 

Ssignal[ri] = Uscale ' Scorrected[n] — rcorrected[n] (Eqn. 9); 

where ascaie is a constant representing the ratio of largest magnitudes of frequency 
components from the water suppressed and non-water suppressed data sets: 

max | Fr{Scorrected[n}} \ 

ascaie = (bqn. 10) 

max | FT{rcorrected[n]} \ 

15 and Fr signifies the Fourier Transform or other similar technique, as will be discussed, 
for computation of frequency from time domain. 

A summation of the independently processed receiver coil data sets is performed to yield 
a single superset of MRS results that takes into account each coil's respective weight, 
20 such that: 

Scombined[k] = ^ Sresults[k] • Wi (Eqn. 1 1 ) 

L 

The absorption spectrum 134, S com bined[k] is then plotted for the combined real 
values of the complex Fourier Transform of the phase-corrected, water-suppressed signal 
with residual water removed. As such, a ID absorption spectrum is displayed on a 
25 monitor or other area interface that is an assimilation of the data sets from each of the 
multiple receiver coils influenced by each coil's weighting. The absorption spectrum 
shows the chemical shift of various metabolites in parts per million (ppm) and can be 
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used for the quantification of metabolites in vivo. The resulting ID spectral plot includes 
peaks from which a clinical diagnosis may be made without a separate interpretation of 
the results of each receiver coil. 

5 For ID plots representing the absorption spectrum, the peak area or the peak 

height representing a particular frequency component, can be calculated and compared to 
other peaks. For example, referring to Fig. 3, a ID plot 140 representing a combined 
absorption spectrum for multiple receiver channel data having a number of peaks is 
shown. While such a plot is not in and of itself novel, to produce such a single plot from 

10 multiple receive channels as described herein is novel. When reviewing the ID plot or 
spectrum 140, a medical professional often uses a creatine peak 142 as a reference for in 
vivo MRS studies. In some cases, the mere presence of a particular peak may indicate a 
clinically significant diagnosis. For example, lactate 144 which, in the example, shows 
up as a doublet, should never occur in a normal human brain. Therefore, the presence of 

15 lactate 144 would potentially indicate damage from acute stroke. As a further example, if 
the NAA peak 146 was lower and the creatine peak 142 and the choline peak 148 were 
higher, the plot 140 would be indicative of a patient stricken with Alzheimer's disease. 

Additionally, beyond the Fourier transform, several other nonparametric 
20 techniques can be used for generating results from MRS studies. Two nonparametric 
techniques that can be used for MRS analysis, which employ adaptive filter bank 
approaches, are the ID Capon method and the ID APES method. The filter bank 
approach to power spectral estimation employs a band-pass filter with a fixed bandwidth 
to estimate the power spectral density (PSD) of a given frequency component. The basic 
25 periodogram, as described in Eqn. 11, is a filter bank approach based on the standard 
Fourier transform. 

N 
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Capon and APES are filter bank approaches that improve the estimate of the PSD 
by creating one data-dependent band-pass filter for each spectral point being estimated. 
The band-pass filter is created by a least-squares minimization process, which attempts to 
minimize the total output power of the filter, yet pass the frequency component of interest 
5 unaltered. A second least-squares process is then implemented to estimate the amplitude 
of the filtered signal. Capon and APES provide more accurate spectral estimates with 
lower sidelobes and narrower spectral peaks than the Fourier Transform based 
periodogram techniques. 

10 However, ID Capon and ID APSES methods fail to take into account, at least in 

a direct way, the damping associated with each signal component in Eqn. 1 . Thus, since 
the detected amplitude will be altered by the damping, it is difficult to accurately estimate 
the amplitudes of the various signal components. To alleviate this problem, 2D Capon 
and 2D APES techniques can be used to obtain a 2D absorption spectrum plot. 

15 

One benefit of a high-resolution 2D spectral plot derived by Capon or APES is 
that peaks that are close in frequency may be separated if the associated damping factors 
are different from one another. Another benefit of the 2D spectral plots is that the 
damping information provided by Capon or APES provides an effective means of 

20 estimating T 2 *, the effective decay rate of transverse magnetization, for several 
metabolites of interest. Furthermore, unlike ID methods, the true amplitudes are not 
masked by damping effects. Empirical results have shown that Capon analysis provides a 
more accurate estimate of frequency and damping since it provides fine resolution to 
resolve closely spaced peaks, while APES provides a more accurate estimate of the 

25 amplitude for a given frequency and damping factor. As will be shown, 2D absorption 
plots may also benefit from the current invention through its multiple receiver channel 
weighting and combination techniques 
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For 2D Capon and 2D APES the filter is subject to the constraints of a damped 
sinusoid of interest, represented by: 

Ae /e e ( - < ™ (Eqn. 13) 

This constraint is based upon the extent that the sinusoid of interest passes undistorted 
5 through the filter, while the total output energy of the filter is minimized. Two- 
dimensional Capon and 2D APES provide an estimate of both the amplitude and phase 
for each frequency and damping factor of an input signal. The 2D Capon and 2D APES 
methods provide high-resolution 2D MRS results showing both frequency, co , and 
damping, <r . In the case of a 2D proton spectroscopy, the 2D Capon method and the 2D 
10 APES method, each improves upon the results of ID proton spectroscopy. The results 
have more defined peaks and more precise amplitude. Additionally, a damping factor 
representing the exponential decay of the associated metabolite is also provided and is 
useful for clinical diagnosis. 

15 The 2D Capon and APES methods may be summarized as follows. Given a 

complex signal, x[n] , « = 0, 1, . . . iV +M - 2, defined: 

x[n] = [x[n] x[n+l] — x[n+M-l]f (Eqn. 14), 

which has a covariance matrix: 

^ M (Eqn. 15), 

71=0 

20 where ( ) H denotes a Hermitian (complex conjugate) transpose. The signal x[n] is 
filtered to produce x aQ) [n] at each a andta of interest. The filtered signal x a & [n] is given 
by: 

x„[n] = h H (a 9 a>)x[»l " = 0,1,-^-1 (Eqn. 16), 

where h (a, co) is an M- vector of FIR filter coefficients. The 2D Capon and 2D APES 
25 filters are specified, respectively, by: 

h Cap on(°,«>) = - H , - (Eqn. 17) 

15 
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and 



where: 



5 and 



where: 



and 



W».»>- S^^ff ■ (Eqn. 18), 

5 (<T,<B)g 0 (<7,ft))s(<7,<0) 



s(a,a>) = [\ e ( -° +Uo) ■■■ e ( -° +J<0 * M - l) ] r (Eqn. 19) 



Q 0 (cj,w) = R, --\-X N (o,co)X N H (cj,(o) (Eqn. 20), 



X N (a,co) = Y 4 x[n}e° n e ] '* n (Eqn. 21) 

n=0 



10 ^ (CT)= ___ (Eqn. 22). 

The least-squares estimate of the spectrum S(g,(d) is then determined at each a andca of 
interest from the filtered signal x a a) [n] by the formula 

S(fr 9 a>) = T ^E^,J«K (T V^ (Eqn. 23) 

or 

15 S(a,a>) = — l —h H {cj,G>)X N {cj,G>) (Eqn. 24); 

where the choice of h{a,(o) determines the method used. 



Referring again to Fig. 2, MRS data is again acquired 110 from multiple receive 
channels of an MR system simultaneously. After data acquisition, the raw data is pre- 
20 processed. That is, once MRS data is acquired 110, non-water suppressed reference data 
is collected, and averaging 114 is used to yield an averaged set of reference data, r[w], 
for each receiver coil, where n represents normalized time. 
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The reference data, r[n], is subjected to the DC mixing 116, zero phasing 118, 
linear phase correction 120, and phase spline smoothing 122 to compute 124 a phase 
correction vector, c[n]. Additionally, the reference data sets, r[n], are used to determine 
5 whether data from a particular coil should be used 125 to generate combined results. If it 
is determined that data from a particular coil should be used, the reference data set, r[n], 
is also used to assign weights 125 to each receiver coil. 

The weightings 125 are used to increase the influence of particular data sets on a 
10 combined superset of data. Specifically, those channels with greater signal-to-noise 
ratios are more heavily weighted. Therefore, data received from regions of interest close 
to receive coil elements are appropriately weighted to take advantage of the improved 
signal-to-noise ratios. The signal characteristics from each coil depend on a number of 
factors including the orientation of the coil with respect to the B 0 field, the proximity of 
15 the coil to the volume generating the signal, coil loading, coil-to-coil coupling effects, 
and the permeability and permitivity of the medium through which the radio-frequency 
signal travels prior to being received by the coil elements. 

A weighted least-squares estimator for 2D Capon or 2D APES may also be 
20 implemented. Specifically, an exponential time-weighting factor to adjust the least- 
squares error minimization criterion for both 2D Capon and 2D APES is incorporated in 
such a way that data samples are weighted differently based on how early or late they 
occur in the time sequence. Thus, signal elements are preferentially weighted in time 
sequentially to take advantage of improved signal-to-noise characteristics earlier in the 
25 signal. The exponential time weighting factor used to adjust the least-squares error 
minimization criterion in 2D Capon and 2D APES is a novel technique used for 
processing MRS data and is thus part of the invention. 
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As previously described, water-suppressed data sets are also collected 110 from 
each coil element and the data from each coil element is averaged to obtain a number of 
averaged, water-suppressed data sets, s[n] 112. Data acquired during the reference 
acquisition is used to compute a phase correction vector 124 that is applied to the water- 
5 suppressed data 126. The reference data is also used to remove residual water. The 
phase correction vector, c[n] 128, is applied to the reference data, r[n] 127, and the 
water-suppressed data, s[n] 126. Applying c[n] removes residual water signal from the 
water-suppressed data sets, r[n]. Water subtraction 128, in which a scaled version of the 
reference data sets is subtracted from the water-suppressed data sets, is performed and a 

10 nonparametric technique, implementing an adaptive filter bank, is used to generate the 
results from the phase-corrected, water-suppressed data sets with residual water effects 
removed. Specifically, either a 2D Capon or 2D APES technique is employed. The 
current invention performs pre-processing of raw data prior to 2D Capon or 2D APES by 
applying the previously described method of MRS data pre-processing involving phase 

1 5 correction and residual water removal. 

In implementing 2D Capon, a frequency spectrum co versus damping cr , 
C results [a,co] , is calculated 150 from the time domain data representing the corrected and 

enhanced MRS signal, S^gnaiM, using the non-parametric filter banks according to the 
20 Capon technique describe above. On the other hand, if 2D APES is implemented, a 
frequency spectrum co versus damping a , A rauUs [(T 9 (o] 9 is calculated 152 from the time 

domain data representing the corrected and enhanced MRS signal, S s \ gm \[n], using the 
non-parametric filter banks according to the APES technique described above. 

25 The frequency spectrum co versus damping a results for each receiver coil, 

including each receive coil's respective weighting, as previously described, are then 

combined 132. If a 2D Capon technique is employed the combination results in a 
supserset of data 154, C combined [a,co] , defined by: 
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^combined 

[a,G)] = 2;C >BWto [<7,fi)]-w f (Eqn. 25). 

L 

However, if a 2D APES technique is employed, the combination results in a superset of 
data 156, A combined [ex, ©], defined by: 

5 From the 2D Capon method and 2D APES method, 2D absorption spectrums, 

CcombinedfrM 138 md A combined\? M 140 > displayed, respectively, for data acquired 
over multiple receiver channels. These absorption spectra allow a reviewing medical 
professional to interpret the data received by each receiver coil simultaneously thereby 
increasing diagnosing time and accuracy. 

10 

Referring now to Fig. 4, a combined 2D absorption spectrum 158 is shown. When 
compared with the ID absorption spectrum of Fig. 3, Fig. 4 shows improved peak-to- 
noise ratios vs. the MRS results from Fig. 3. The ratios of peak heights vs. the standard 
deviation of a noise region (near 0.0 ppm) are provided for several peaks of interest. This 

15 distinction becomes more evident when comparing peak-to-noise ratios for 2D Capon or 
2D APES with each peak of a corresponding noise region with the same value of 
damping, a. Again, by compiling a single absorption spectrum, rather than requiring a 
medical professional to review an absorption spectrum from each receiver coil 
independently and individually, the burden on the reviewing medical professional is 

20 lowered and typically makes the review process considerably more efficient. The 
combined 2D absorption spectrum 158 allows a reviewing medical professional to 
interpret the data received by each receive coil simultaneously. 

The invention described above is a technique for multi-channel MRS whereby 
25 MR signals are simultaneously acquired from multiple coils and individually processed to 
generate multiple sets of MRS results. Furthermore, the multiple sets of MRS results are 
combined to form a single superset of MRS results, such that the MRS results are 
displayed as a proton MRS absorption spectrum from the single superset of MRS results. 

19 
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It is contemplated that the above architecture and technique can be embodied in any 
modality of a medical imaging scanner capable of providing data required for MRS. 

It is further contemplated that the above technique be embodied in a MRI system 
5 having a plurality of gradient coils positioned about a bore of a magnet to impress a 
polarizing magnetic field and an RF transceiver system and an RF switch controlled by a 
pulse module to transmit RF signals to an RF coil assembly to acquire MR images. The 
MRI system also includes a computer programmed to acquire MR signals from multiple 
receiver channels and process the MR signals from each channel independently. The 
10 computer then combines the processed MR signals from each coil into a single set of 
MRS results and displays a proton MRS absorption spectrum from the single composite 
set of MRS results. 

It is additionally contemplated that the above technique be embodied in a 
15 computer program including instructions which when executed by a computer cause the 
computer to process MR data from the multiple coils independently. The computer is 
also caused to combine the processed MR data into a single composite set of MRS results 
and display the set as a single MRS absorption spectrum. 

20 The present invention has been described in terms of the preferred embodiment, 

and it is recognized that equivalents, alternatives, and modifications, aside from those 
expressly stated, are possible and within the scope of the appending claims. 
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